A=[1 1 1 1 1;-1 1 -1 1 -1;2 4 8 16 32;-2 4 -8 16 -32;3 9 27 81 243];
for i=2:5
A(:,i)=A(:,i)./factorial(i);
end
format rat
C=inv(A)
M=[1 0 0 0 0;0 -1 0 0 0;-1 0 1 0 0;0 1 0 -1 0;0 0 -1 0 1];
format rat
T5=C*inv(M)
